A ndreev- Tunneling, Coulomb Blockade, and 
Resonant Transport of Non-Local Spin-Entangled Electrons 
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We propose and analyze a spin-entangler for electrons based on an s-wave superconductor coupled 
to two quantum dots each of which is tunnel-coupled to normal Fermi leads. We show that in the 
presence of a voltage bias and in the Coulomb blockade regime two correlated electrons provided by 
the Andreev process can coherently tunnel from the superconductor via different dots into different 
leads. The spin-singlet coming from the Cooper pair remains preserved in this process, and the 
setup provides a source of mobile and nonlocal spin-entangled electrons. The transport current is 
calculated and shown to be dominated by a two-particle Breit-Wigner resonance which allows the 
injection of two spin-entangled electrons into different leads at exactly the same orbital energy, which 
is a crucial requirement for the detection of spin entanglement via noise measurements. The coherent 
tunneling of both electrons into the same lead is suppressed by the on-site Coulomb repulsion and/or 
the superconducting gap, while the tunneling into different leads is suppressed through the initial 
separation of the tunneling electrons. In the regime of interest the particle-hole excitations of the 
leads are shown to be negligible. The Aharonov-Bohm oscillations in the current are shown to 
contain single- and two-electron periods with amplitudes that both vanish with increasing Coulomb 
repulsion albeit differently fast. 



I. INTRODUCTION 

The creation of nonlocal pairwise-entangled quantum 
states, so-called Einstein-Podolsky-Rosen (EPRA pairsEl, 
is essential for secure quantum cammunicationa, dense 
coding and quantum teleportatiorjj, or more fundamen- 
tal, for testing the violation of Bell's inequality!! Such 
tests already exist for photons but not yet for massive 
particles such as electrons since it is difficult to pro- 
duce and to detect entangled electrons. However, there 
is strong experimental evidence that electron spins in a 
semiconductor environment show unusually long dcphas- 
ing times approaching microseconds and that they can 
be transpartedj-ohase coherently over distances exceed- 
ing lOO/incmoiJ. This makes spins of electrons in semi- 
conductors promisinftjcandidates for carriers of quantum 
information (qubits)Ei||13. In particular, we have recently 
proposed a setupEj consisting of a spin-entangler and a 
beam splitter where the spin-entanglement is detectable 
via electronic transport properties. We have shown that 
the current-current correlations (noise) is enhanced if the 
entangled electrons are spin-singlets leading to bunching 
behavior, whereas the noise is suppressed for spin-triplets 
leading tp-antibunching behavior. 

In RefO we assumed the existence of an entangler, i.e. 
a device that generates spin-singlets which are made out 
of two electrons which reside in different but degenerate 
orbital states, and we focussed on the question of how to 
detect spin-entangled electrons via transport and noise 
measurements. Here, we address the problem of how to 
implement pich an entangler in a solid state device. We 
have foundES that for such noise measurements, which 
are based on two-particle interference effects, it is abso- 



lutely crucial that both electrons, coming from different 
leads, possess the same orbital energy. If the orbital en- 
ergies of the two entangled electrons are different, the 
electrons cannot interfere with each other, and thus spin 
correlations would not be observable in the noise£3. 




FIG. 1. The entangler setup: Two spin-entangled electrons 
forming a Cooper pair can tunnel with amplitude Tsd from 
points ri and T2 of the superconductor, SC, to two dots, Di 
and L>2, by means of Andreev tunneling. The dots are tun- 
nel-coupled to normal leads Li and L2, with tunneling ampli- 
tude Tdl- The superconductor and leads are kept at chemical 
potentials fis and fii, resp. 

In the following we propose a setup which involves 
a superconductor coupled to two quantum dots which 
themselves are coupled to normal leads, see Fig. 1. We 
show that such a setup acts as an entangler which meets 
all the requirements needed for a successful detection of 
spin .Entanglement via noise measurements. In previous 
worlo we showed that in equilibrium the spin correla- 



tions of an s-wave superconductor induce a spin-singlet 
state between two electrons each of which resides on a 
separate quantum dot which both are weakly coupled to 
the same superconductor (but not among themselves). 
This non-local spin-entanglement leads then to obsept- 
able effects in a generalized Josephson junction setupEJ. 
In the present work we consider a non- equilibrium situ- 
ation where an applied voltage bias drives a stationary 
current of pairwise spin-entangled electrons from the su- 
perconductor through the quantum dots into the leads, 
see Figs. 1 and 2. 



II. QUALITATIVE DESCRIPTION OF THE 
ANDREEV ENTANGLER 

We begin with a qualitative description of the entan- 
gler and its principal mechanism based on Andreev pro- 
cesses and Coulomb blockade effects. In subsequent sec- 
tions we introduce then the Hamiltonian and calculate 
the stationary current in detail. We consider an s-wave 
superconductor which acts as a natural source of spin- 
entangled electrons, since the electrons form Cooper pairs 
with singlet spin-wavefunctionsla. The superconductor, 
which is held at the chemical potential fj,s, is weakly cou- 
pled by tunnel barriers to two separate quantum dots D\ 
and -D2 which themselves are weakly coupled to Fermi 
liquid leads L\ and L2, resp., both held at the same 
chemical potential /ii — fj.2- The corresponding tunneling 
amplitudes between superconductor and dots, and dots- 
leads, are denoted by T$d and Tdl, resp. (for simplicity 
we assume them to be equal for both dots and leads). 

In general, the tunnel- coupling of a superconductor to 
a normal region allows for coherent transport of two elec- 
trons of opposite spins due to Andreev .tunnclingjla, while 
single-electron tunneling is suppressed!^. In the present 
setup, we envision a situation where the two electrons 
are forced to tunnel coherently into different leads rather 
than both into the same lead. This situation can be 
enforced in the presence of two intermediate quantum 
dots which are assumed to be in the Coulomb blockade 
rcgimco so that the state with the two electrons being on 
the same quantum dot is strongly suppressed, and thus 
the electrons will preferably tunnel into separate dots and 
subsequently into separate leads (this will be quantified 
in the following). 

By applying a bias voltage A/i = us — \Ai > transport 
of entangled electrons occurs from the superconductor 
via the dots to the leads. The chemical potentials ei 
and £2 of the_jquantum dots can be tuned by external 
gate voltagesEa such that the coherent tunneling of two 
electrons into different leads is at resonance, described 
by a two-particle Breit-Wigner resonance peaked at t\ = 
e 2 = fJ-s- In contrast, we will see that the current for the 
coherent tunneling of two electrons into the same lead 
is suppressed by the on-site Coulomb U repulsion of a 
quantum dot and/or by the superconducting gap A. 



Next, we introduce the relevant parameters describing 
the proposed device and specify their regime of inter- 
est. First we note that to avoid unwanted correlations 
with electrons already on the quantum dots one could 
work in the cotunncling regimeEJ where the number of 
electrons on the dots are fixed and the resonant levels 
e;, I = 1,2 cannot be occupied. However, we prefer to 
work at the resonances e/ — /^Sj since then the total cur- 
rent and the desired suppression of tunneling into the 
same lead is maximized in this regime. Also, the de- 
sired injection of the two electrons into different leads 
but at the same orbital energy is then achieved. It turns 
out to be most convenient to work in the regime where 
the dot levels ei have vanishing occupation probability. 
For this purpose we require that the dot-lead coupling 
is much stronger than the superconductor-dot coupling, 
i.e. \Tsd\ < \Tdl\: so that electrons which enter the dots 
from the superconductor will leave the quantum dots to 
the leads much faster than new electrons can be provided 
from the superconductor. In addition, a stationary occu- 
pation due to the coupling to the leads is exponentially 
small if A/i > koT, T being the temperature and ko 
the Boltzmann constant. Thus in this asymmetric bar- 
rier case, the resonant dot levels e/ can be occupied only 
during a virtual process. 

Next, the quantum dots in the ground state are allowed 
to contain an arbitrary but even number of electrons, 
No = even, with total spin zero (i.e. antifcrromagnetic 
filling of the dots). An odd number No must be excluded 
since a simple spin-flip on the quantum dot would be pos- 
sible in the transport process and as a result the desired 
entanglement would be lost. Further, we have to make 
sure that also spin flip processes of the following kind 
are excluded. Consider an electron that tunnels from 
the superconductor into a given dot. Now, it is possi- 
ble in principle (e.g. in a sequential tunneling processed) 
that another electron with the opposite spin leaves the 
dot and tunnels into the lead, and, again, the desired 
entanglement would be lost. However, such spin flip pro- 
cesses will be excluded if the energy level spacing of the 
quantum dots, Se, (assumed to be similar for both dots) 
exceeds both, temperature fc^T and bias voltage A/i. A 
serious source of entanglement-loss is given by electron 
hole-pair excitations out of the Fermi sea of the leads 
during the resonant tunneling events. However, we show 
in the following that such many-particle contributions 
are suppressed if the resonance width 7; = 2tti>i\T£)l\ 2 
is smaller than A/x (for e; ~ /is) 7 where vi is the den- 
sity of states (DOS) per spin of the leads at the chemical 
potential \ii. 

Finally, an additional energy scale that enters the con- 
sideration is the superconducting gap energy A, which is 
half the minimum energy it costs to break up a Cooper 
pair into two quasiparticles. This gap energy also char- 
acterizes the time delay between the subsequent coherent 
Andreev tunneling events of the two electrons of a Cooper 
pair. In order to exclude single-electron tunneling where 
the creation of a quasiparticle in the superconductor is a 



final excited state we require that A > A/j,, ksT. 

To summarize, the regime of interest in this work is 
then given by 

A, U, Se > A^ > 71, k B T, and ji > j s ■ (2.1) 

Some inequalities will become clear when we discuss the 
various processes in detail below. As regards possible ex- 
perimental implementations of the proposed setup and 
its parameter regime, we would like to mention that, 
typically, quantum dots are made out of semiconducting 
heterostructures, which satisfy above inequalitiesc-1 Fur- 
thermore, in recent experiments, it has been shown that 
the fabrication of hybrid structures with semiconductor 
and superconductor being tunnel- coupled is possibleEj'Ej. 
Other candidate materials are e.g. carbon nanotubes 
which also show Coulomb blockade behawor with U and 
Se being in the regime of interest herdlJ. The present 
work might provide further motivation to implement the 
structures proposed here. 

Our goal in the following is to calculate the stationary 
charge current of pairwise spin-entangled electrons for 
two competing transport channels, first for the desired 
transport of two entangled electrons into different leads 
(ii) and second for the unwanted transport of both elec- 
trons into the same lead (12)- We compare then the two 
competing processes and show how their ratio, I1/I2, de- 
pends on the various system parameters and how it can 
be made large. An important finding is that when tun- 
neling of two electrons into different leads occurs, the cur- 
rent is suppressed due to the fact that tunneling into the 
dots will typically take place from different points ri and 
T2 on the superconductor (see Fig. 1) due to the spatial 
separation of the dots D\ and £>2. We show that the dis- 
tance of separation Sr — \r\ — r2 1 leads to an exponential 
s uppr ession of the current via different dots if Sr > £ (see 
( 5.13 )), where £ is the coherence length of a Cooper pair. 
In the relevant regime, Sr < £, however, the suppression 
is only polynomial and ex 1/ \kp5r) 2 , with Izf being the 
Fermi wavevector in the superconductor. On the other 
hand, tunneling via the same dot implies Sr = 0, but suf- 
fers a suppression due to U and/or A. The suppression 
of this current is given by the small parameter {ji/U) 2 
in the case U < A, or by (ji/A) 2 , if U > A as will be 
derived in the following. Thus, to maximize the efficiency 
of the entangler, we also require kpSr < A, U. 

Finally, we will discuss the effect of a magnetic flux on 
the entangled current in an Aharonov-Bohm loop, and 
we will see that this current contains both, single- and 
two-particle Aharonov-Bohm periods whose amplitudes 
have different parameter dependences. 



H = H s +J2H Dl +J2 H Li, 1=1,2. (3.1) 

1 1 

Here, the superconductor is described by the BCS- 
Hamiltonianli3 Hg — J2 ^k7k CT 7kcr, where a =1,1, and 

k,o- 

the quasiparticle operators 7k<j describe excitations out 
of the BCS-groundstate |0)s defined by 7kir|0),s = 0. 
They are related to the electron annihilation and cre- 
iCko- and c ko . through the Bogoliubov 



ation operator* 
transformation^ 



CkT = Wk7kT + v kJ- k i 
c-ki = "k7-ki - «k7kT 



(3.2) 



where w k = (1/V2)(1 + £ k /£ k ) 1/2 and v w = (l/yS)(l - 
Ck/^k) 1 ^ 2 are the usual BCS coherence factorsEa, and 
£k = ^k ~ Ms is the normal state single-electron energy 
counted from the Fermi level fxg, and Ey. = ^Ck + ^ 2 
is the quasiparticle energy. We choose energies such that 
Us = 0. Both dots are represented as one localized 
(spin-degenerate) level with energy e; and is modeled 
by an Anderson- type Hamiltonian Hm — e; ^ w a di a + 

a 

Uni^nii, I = 1, 2. The resonant dot level ei can be tuned 
by the gate voltage. Other levels of the dots do not par- 
ticipate in transport if Se > A/x > fcgT, where A/i = 
— ijli, and ni is the chemical potential of lead I = 1,2, and 
Se is the single particle energy level spacing of the dots. 
The leads I = 1,2 are assumed to be non-interacting 
(normal) Fermi liquids, B. u = J2kcr e ^ a lka a i^- Tunnel- 
ing from the dot I to the lead / or to the point r; in the 
superconductor is described by the tunnel Hamiltonian 
Ht = H S d + H DL with 



H SD =Y,TsDd\Mv l )+^ 

la 

Hdl = ^2 T DLal ka di a + h.c. . 



(3.3) 

(3.4) 



Iko 



Here, ^v(rj) annihilates an electron with spin a at site 
Yi, and d] creates it again (with the same spin) at dot 
I with amplitude T$d- ^cr(rj) is related to Cko- by the 
Fourier transform ip a (ri) — ^2e lkri Cka- Tunneling from 

k 

the dot to the state k in the lead is described by the tun- 
nel amplitude T^l- We assume that the k-dependence 
of Tdl can be safely neglected. 



IV. STATIONARY CURRENT AND T-MATRIX 



III. HAMILTONIAN OF THE ANDREEV 
ENTANGLER 

We use a tunneling Hamiltonian description of the sys- 
tem, H — Hq + Ht, where 



The stationary current of two electrons passing from 
the superconductor via virtual dot states to the leads is 
given by 



I = 2eJ2WfiPi, 



(4.1) 



where Wf% is the transition rate from the superconductor 
to the leads. We calculate this transition rate in terms 
of a T-matrix approached, 



W fl = 2Tr\(f\T(e i )\ l )\ 2 5(e f ~e t ). 



(4.2) 



Here, T(ei) = Ht e . + l H {£i — H ), is the on-shell trans- 
mission or T-matrix, with 77 being a small positive real 
number which we take to zero at the end of the calcu- 
lation. Finally, pi is the stationary occupation proba- 
bility for the entire system to be in the state \i). The 
T-matrix T(si) can be written as a power series in the 
tunnel Hamiltonian Ht, 



T(s t ) = H T + H T J2( 



1 



\ £i + if] - H a 



HtT , (4.3) 



where the initial energy is e% — 2p$ = 0. We work 
in the regime defined in Eq. (2.1), i.e. 7; > 75, and 
A,U,8e > A/i > 7;,/cbT, and around the resonance 
e L ~ us. Further, 75 = 2ttv s \Tsd\ 2 and 7/ = 2iti>i\T dl \ 2 
define the tunneling rates between superconductor and 
dots, and between dots and leads, respectively, with vs 
and v\ being the DOS per spin at the chemical poten- 
tials [is and fii , respectively. We will show that the total 
effective tunneling rate from the superconductor to the 
leads is give n b y 75/7/ due to the Andreev process. In 
the regime (2.1) the entire tunneling process becomes a 
two-particle problem where the many-particle effect of 
the reservoirs (leads) can be safely neglected and the co- 
herence of an initially entangled Cooper pair (spin sin- 
glet) is maintained during the transport into the leads 
as we shall show below. Since the superconducting gap 
satisfies A > A/i, ksT, the superconductor co ntai ns no 
quasiparticle initially. Further, in the regime ( |2.l| ), the 
resonant dot levels t\ are mostly empty, since in the as- 
sumed asymmetric case, \Tdl\ > \Tsd\ (or 7/ > 7s), the 
electron leaves the dot to the lead much faster than it can 
be replaced by another electron from the superconductor. 
In addition, we can neglect any stationary occupation of 
the dots induced by the coupling of the dots to the leads. 
Indeed, in the stationary limit and for given bias A/i 
this occupation probability is determined by the grand 
canonical distribution function oc exp(— A/i/fe^T) <C 1, 
and thus pi ~ for any initial state where the resonant 
dot level is occupied. In this regime, the initial state \i) 
becomes \i) — \0)s\0)d\Hi)i, where \0)s is the quasiparti- 
cle vacuum for the superconductor, \0}d means that both 
dot levels e; are unoccupied, and \pi)i defines the occupa- 
tion of the leads which are filled with electrons up to the 
chemical potential p,\. We remark that in our regime of 
interest no Kondo effects appear which could destroy the 
spin entanglement, since our dots contain each an even 
number of electrons in the stationary limit. 



V. CURRENT DUE TO TUNNELING INTO 
DIFFERENT LEADS 

We now calculate the current for simultaneous coher- 
ent transport of two electrons into different leads. The 
final state for two electrons, one of them being in lead 1 
the other in lead 2, can be classified according to their to- 
tal spin S. This spin can be either a singlet (in standard 
notation) \S) = (| U) - | IT))/V2 with S = 0, or a triplet 
with S — 1. Since the total spin is conserved, [S 2 , H] = 0, 
the singlet state of the initial Cooper pair will be con- 
served in the transport process and the final state must 
also be a singlet. That this is so can also be seen ex- 
plicitly when we allow for the possibility that the final 
state could be the S z = triplet |* ) = (| U) + | IT)) A/2- 
[The triplets \t+) = | ft) and |i_) = | |j) can be ex- 
cluded right away since the tunnel Hamiltonian Ht con- 
serves the spin-component a and an Andreev process in- 
volves tunneling of two electrons with different spin a.] 
Therefore we consider final two-particle states of the form 



and 



I/} = (W2)[<4pT a 2q.l ±a lpi a 2qT]N)' Where the 

+ signs belong to the singlet \S) and triplet |t ), resp. 
Note that this singlet/triplet state is formed out of two 
electrons, one being in the p state in lead 1 and with 
energy e p , while the other one is in the q state in lead 
2 with energy e q . Thus, the two electrons are entangled 
in spin space while separated in orbital space, thereby 
providing a non-local EPR pair. The tunnel process to 
different leads appears in the following order. A Cooper 
pair breaks up, where one electron with spin a tunnels to 
one of the dots (with empty level e;) from the point of the 
superconductor nearest to this dot. This is a virtual state 
with energy deficit E^ > A. Since A > 7;, the second 
electron from the Cooper pair with spin — a tunnels to 
the other empty dot-level before the electron with spin a 
escapes to the lead. Therefore, both electrons tunnel al- 
most simultaneously to the dots (within the uncertainty 
time h/A). Since we work at th e resonance e; ~ fxs = 
the energy denominators in (4.3) show divergences ex 1/r/ 
indicating that tunneling between the dots and the leads 
is resonant and we have to treat tunneling to all orders 
in Hdl in (4.3), eventually giving a finite result in which 
77 will be replaced by 7z/2. Tunneling back to the su- 
perconductor is unlikely since \Tsd\ < \Tdl\- We can 
therefore write the transition amplitude between initial 
and final state as 

(f\ T o\i) = ^=(a 2q xai pT T 4 T 4i)([ d 2|dit ± ^21^11]^ ) , 

(5.1) 

where To = T(ei = 0), and the partial T-matrices T and 
T are given by 



T = —H S d- ttH S d 

ir) - H 17] - H 



(5.2) 



and 



T =H 



n=0 



1 



ii] - H a 



-H 



DL 



2ra+l 



(5.3) 



In (5.1) we used that the matrix element containing T is 
invariant under spin exchange t < ~ > l, and the abbreviation 
(...) stands for (i|...|i). The part containing T describes 
the Andreev process, while the part containing T is the 
resonant dot <-> lead tunneling. 

We first consider the Andreev process. We insert a 
complete set of single-quasiparticle (virtual) states, i.e., 
1 = Y J ,v„li„d)_ rr \i)(i\di- a '^ a , between the two H S d 
in (5.2) and use that the resulting energy denominator 
|z?; — i^k — e;| w l-Ekl, since we work close to the resonance 
ei ~ and £k > A. The triplet contribution vanishes 
since Uk«k = w_k«-k for s-wave superconductor. For the 
final state being a singlet we then get 



((^2|diT -d 2 -\d n )T ) 



AT 2 
* 1 SD 



e\ + e 2 - if] 



E ~»kfk 
. E v 



cos (k • 8t) , (5.4) 



where St — ri — r 2 denotes the distance vector between 
the points on the superconductor from which electron 1 
and 2 tunnel into the dots. To evaluate the sum over k 
we use Mk^k = A/(2i?k), linearize the spectrum around 
the Fermi level with Fermi wavevector fop, and obtain 
finally for the Andreev contribution 



(0*2|diT -d 2 -\d n )T 



2-KvsTg D sm(kpSr) 
ei + £2 — if} UfSt 



(5.5) 



A. Dominant contribution of resonant tunneling to 
different leads 



Now we calculate the matrix element in ( |5.1[ ) contain- 
ing T where tunneling has to be treated to all orders 
in Ht- To simplify the notation we suppress spin in- 
dices and introduce a ket notation |12), where 1 stands 
for quantum numbers of the electron on dot 1/lead 1 
and similar for 2. E.g. \pq) stands for a\ a\ a \i), 
where p is from lead 1 and q from lead 2; or \pD) 
stands for o\ d\ _ a \i), etc. We concentrate first on 
the resummation of the following dot «-v lead transitions 
\DD) -> \LD) -> \DD) or \DD) -> \DL) -> \DD). 
In this sequence, \DD) is the state with electron on 
dot 1 and the other one on dot 2, and \LD) defines a 
state where one electron is in lead 1 and the other one 
on dot 2. We exclude processes of the kind \DD) — > 
\LD) -v \LL) -> \LD) -v \DD) or \DD) -> \LD) -> 
\LL) — > \DL) — * \DD), where both electrons are virtu- 
ally simultaneously in the leads as well as the creation of 
electron-hole excitations out of the Fermi sea. We show 
in App. A and B that such contributions are suppressed 



in the regime (2.1) considered here by the small param- 
eter 7;/A/x. The dominant contribution is then given by 

(pq\T'\DD) 

OO 

(pq\H DlLl \Dq)(Dq\ £ (-l^Ho^f^Dq) 

n—0 

^{Dq\ J ^ w H D2L2 \DD) 

OO 

+ {pq\H D2L2 \pD)( P D\ £ (-l^-Hn^^lpD) 

n—0 

y-(pD\ T ^ w H DlLl \DD) 

oo 1 

x{DD\ J2(-- jrH DL f m \DD). (5.6) 

Since the sums for the transition \DD) — ► \DD) via the 
sequences \DD) -> \LD) -> |DD) and |DD) -> |DL) -> 
l-D-D) are independent, we can write all summations in 
( |5.6[ ) as geometric series which can be resummed explic- 
itly. We begin with the two-particle process for which we 
find 



(DD\jr,(^—H DL f m \DD) 



1-(DD\(^H DL )*\DD) 



where 
(DD\ 



1 



v ir] - H 



-H 



DL) 



\DD) = 



ir] - ei - e 2 



(5.7) 



(5.8) 



with E being the self-energy, E = |l£> L | 2 J2ik( ir ) - e i - 
£k) _1 - In the presence of a Fermi sea in the leads, we 
introduce a cut-off in the sum in E at the Fermi level 
£k ~ — A/i and at the edge of the conduction band, e c . 
Then we obtain E = ReE — i-f/2, where 7 = 71+72, 
and the logarithmic renormalization of the energy level 
is small, i.e. ReE ~ 7; ln(e c /A/i) <C A/i and will be 
neglected. Finally, we arrive at the following expression 



(DD\J2(^—H DL f™\DD)= * + *-** 



— ' irj - H 

rn—0 ' 



ei + e 2 — 27/2 



(5.9) 



Simi lar results hold for the one-particle resummations in 



1 



£2 + £ P - *?? 



fr* «??--Ho £2+£ p -«72/2 



(Dq\J2( 



1 



»=0 



27/ - Ho 



H DlLl Y n \Dq) 



(5.10) 
£1 + e q - ir\ 

£l + £q - *7l/2 

(5.11) 



Inserting the preceding results back into Eq. ( p]q ) we 
obtain 



(pq\T \DD) 



- t dl( £ i + £ 2 - ir]) 



(ei + e q - »7i/ 2 )( e 2 + e P - i72/2) ' 

(5.12) 



Thus, we see that the resummations cancel all diver- 
gences like the (ei + £2 — iff) denominator appearing in 
( |5.4D and (|5.5|), and that, as expected, the resummation 
of divergent terms leads effectively to the replacement 
irj — ► iji/2 so that the limit e; — > is well-behaved. 
It is interesting to note that the two-particle resonance 



(ei + €2 — J7/2) occurring in (5.9) has canceled out 



in ( 5.12 ), and we finally obtain a product of two inde- 
pendent single-particle Breit-Wigner resonances. Still, 
we will just see that the two-particle correlation is rein- 
troduced wh en we insert (5.12) into the expression for 
the current (4.1) due to the integrations over p, q, and 
the fact that the main contribution comes from the reso- 



nances. Indeed, making use of Eqs. (4.1,4.2), and energy 



conservation s 



0, 



-e q , and of Eqs. 



(5.5) and ( 5.12| ), we finally obtain for the current (de- 
noted by If) where each of the two entangled electrons 
tunnels into a different lead 



h 



nh 



(ei+e 2 ) 2 +7 2 /4 



sm(kpSr) 



kp5r 



■ 2dr 
ex P{ t\ , 



< 



(5.13) 



where, again, 7 = 71 + 72- We note that Eq. ( 5.13| ) also 
holds for the case with 71 ^ 72. The current becomes 
exponentially suppressed with increasing distance 5r be- 
tween the tunneling points on the superconductor, the 
scale given by the Cooper pair coherence length £. This 
does not pose severe restrictions for conventional s-wave 
material with £ typically being on the order of microme- 
ters. More severe is the restriction that kpSr should not 
be too large compared to unity, especially if k~f of the 
superconductor assumes a typical value on the order of 
a few Angstroms. Still, since the suppression in kpSr is 
only power-law like there is a sufficiently large regime on 
the nanometer scale for Sr where the current I\ ca n as- 
sume a finite measurable value. The current (5.13) has 



again a Breit-Wigner resonance form which assumes it 
maximum value when e\ = —62, 



h = 



4e7| 
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(5.14) 



This resonance at t\ = —€2 clearly shows that the cur- 
rent is a correlated two-particle effect (even apart from 
any spin correlation) as we should expect from the An- 
dreev process involving the coherent tunneling of two 
electrons. Together with the sing le-particle resonances 
discussed above (see after Eq. ( |5.12| )) we thus see that the 
current is carried by correlated pairs of electrons whose 
energies satisfy |e p | = |e q | < 7, if £1 = £ 2 = 0. 



A particularly interesting case occurs when the ener- 
gies of the dots, ei and £2, are both tuned to zero, i.e. 
ei = e 2 = us = 0. We stress that in this case the elec- 
tron in lead 1 and its spin-entangled partner in lead 2 
have exactly the same orbital energy. We have shown 
previouslyEj that this degeneracy of orbital energies is a 
crucial requirement for noise measurements in which the 
singlets manifest themselves in form of enhanced noise in 
the current (bunching), whereas uncorrelated electrons, 
or, more generally, electrons in a triplet state, lead to a 
suppression of noise (antibunching) . 

We remark again that the current I± is carried by elec- 
trons which are entangled in spin space and spatially 
separated in orbital space. In other words, the station- 
ary current I\ is a current of non-local spin-based EPR 
pairs. Finally, we note that due to the singlet character 
of the EPR pair we do not know whether the electron 
in, say, lead 1 carries an up or a down spin, this can 
be revealed only by a spin- measurement. Of course, any 
measurement of the spin of one (or both) electrons will 
immediately destroy the singlet state and thus the entan- 
glement. Such a spin measurement (spin read-out) can be 
performed e.g|— by making use of the spin filtering effect of 
quantum dotd 2 .! The singlet state will also be destroyed 
by spin-dependent scattering (but not-by Coulomb ex- 
change interaction in the Fermi seatj). However, it 
is known experimentally that electron spins in a semi- 
conductor environment show unusually long dephasing 
times approaching microseconds and can be trar. 
phase coherently over distances exceeding lOO/ir 
This distance is sufficiently long for experiments per- 
formed typicalUi-|On the length scale of quantum confined 
nanostructureso. 



B. Negligible tunnel contributions 

We turn now to a discussion of various tunnel pro- 
cesses which we have not taken into account so far and 
show that they are negligibly small compared to the 
ones we have retained. As we mentioned above we ex- 
clude virtual states where both electrons are simultane- 
ously in the leads. This is justified in the regime (2.1) 
considered here. To show this we consider the process 
\DD) — > \DD). This transition occurs either in a tran- 
sition sequence of the type \DD) — > \LD) — > \DD), as 
considered above, leading to the amplitude Adl = — *7i 
(see Eq. (Al) in App. A), or in a sequence of the type 
\DD) -v \LD) -> \LL) -> \DL) -> \DD), where both 
electrons are simultaneously in the leads (|LL)-state), 



1l 



^ + ln (^J 



leading to the amplitude All 

(see Eq. ( |A3| ) in App. A). However, this amplitude All 
is suppressed by a factor 7//A/1 < 1 compared to Adl- 
Further, a process where we create an electron-hole pair 
out of the Fermi sea of the leads could, in principle, de- 
stroy the spin-correlation of the entangled electron pair 
when an electron with the "wrong" spin (coming from the 



Fermi sea) hops on the dot. But such contributions cost 
additional energy of at least A^i, and again such particle- 
hole processes are suppressed by a factor (•yi/Afi) 2 as we 
show in detail in Appendix B. 



VI. TUNNELING VIA THE SAME DOT 

The two electrons of a Cooper pair can also tunnel via 
the same dot into the same lead. In this section we calcu- 
late the current induced by this process. We show that we 
obtain a suppression of such processes by a factor {n/U) 2 
and/or (-yi/A) 2 compared to the process discussed in the 
preceding section. However, in contrast to the previous 
case, we do not get a suppression resulting from the spa- 
tial separation of the Cooper pair on the superconductor, 
since here the two electrons tunnel from the same point 
either from ri or r2 (see Fig. 2). As before, a tunnel pro- 
cess starts by breaking up a Cooper pair followed by an 
Andreev process with two possible sequences, see Fig. 2. 
(I) In a first step, one electron tunnels from the super- 
conductor to, say, dot 1, and in a second step the second 
electron also tunnels to dot 1. There are now two elec- 
trons on the same dot which costs additional Coulomb 
repulsion energy U, thus this virtual state is suppressed 
by 1/U. Finally, the two electrons leave dot 1 and tun- 
nel into lead 1. (II) There is an alternative competing 
process (which avoids the double occupancy). Here, one 
electron tunnels to, say, dot 1, and then the same elec- 
tron tunnels further into lead 1, leaving an excitation on 
the superconductor which costs additional gap energy A 
(instead of U) , before finally the second electron tunnels 
from the superconductor via dot 1 into lead 1. 

We first concentrate on the tunneling process (II), and 
note that the leading contribution comes from the pro- 
cesses where both electrons have left the superconductor 
so that the system has no energy deficit anymore. We 
still have to resum the tunnel processes from the dot to 
the lead to all orders in the tunnel Hamiltonian Hdl- 
In what follows we suppress the label I = 1,2 since the 
setup is assumed to be symmetric and tunneling into ei- 
ther lead 1 or lead 2 gives the same result. The transition 
amplitude (f\To\i) including only leading terms is 

{f\T \i)=Y,{f\ H DL\Dp"a) 

p"<7 
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x(ZVV| £(-__# DL f' W ' CT ) 



:(Dp"a\ 



K irj - H 

-Hsd . 1 TT H DL . l TT H SD \i), (6.1) 
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where again |/) = (l/v^a^a^ ± a^a^^i), with 
± denoting the triplet (+) and singlet (— ), resp., and 
the intermediate state \Dp"a) = d_ 



a u p"a 



i). The index 



a appearing together with the tunnel Hamiltonians in 
(S.l) determines the spin of the electron that tunnels. 





FIG. 2. Two competing virtual processes are shown when 
the two electrons tunnel via the same dot: (a) Andreev pro- 
cess leading to a double occupancy of the dot with virtual en- 
ergy 1/U , and (b) the process which differs by the sequence of 
tunneling leading to an additional virtual energy 1/A instead 
of 1/U. 

There are some remarks in order regarding Eq. ( p.l| ). 
The electron which tunnels to the state |p"er) has not 
to be resummed further since this would lead either to a 
double occupancy of the dot which is suppressed by l/U, 
or to the state with two electrons simultaneously in the 
lead with a virtual summation over the state p". But 
we already showed that the latter p roces s is suppressed 
by 7//A^i. Making then use of Eq. ( |5.1l| ), we obtain for 
the first factor in ( |6.l[ ) 
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where again in ( |6.2| ) and (6.3) the upper sign belongs to 
the triplet and the lower sign to the singlet. For the third 



line of (6.1) containing the superconductor-dot transi- 
tions we obtain 
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Combining the results (p.2|)-( |6.4| ) we obtain for the am- 
plitude (O) 



(f\To\i) 
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for the final state |/) being a singlet, whereas we get 
again zero for the triplet. 

Next we consider the process (I) where the tunneling 
involves a double occupancy of the dot (see left panel 
in Fig. 2). In this case the transition amplitude can be 
written as 



(/|T |i) = J2(f\ H DL\Dp"a) 
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As before, the transition amplitude (f\To\i) is only 
nonzero for the final lead state |/) being a singlet state. 
Repeating a similar ca lcula tion as before we find that the 
amplitude is given by ( |6.5|) but with A bei ng r eplaced b y 
U/ir. We note that the two amplitudes ( |6.5D and ( |6.6| ) 
have the same initial and same final states. Thus, to ob- 
tain the total current due to processes (I) and (II ) we 
need to add these two amplitudes. Then, using Eq. (4.1) 
we find for the total current I2 in case of tunneling of two 
electrons into the same lead, 
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We see that the effect of the quantum dots consists in 
the suppression factor (j/£) 2 for tunneling into the same 
lead. We remark that in contrast to the previous case 
(tunneling into different leads) the current does not have 
a resonant behavior since the virtual dot states are no 
longer at resonance due the energy costs U or A in the 
tu nnelin g process. Our final goal is to compare I\ given 
in (5.14) with Ii- Thus, forming the ratio of the currents 
of the two competing processes, we obtain 
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From this ratio we see that the desired regime with I\ 
dominating I2 is obtained when f/7 > kp5r, and Sr < £. 



We would like to emphasize that the relative suppression 
of I2 (as well as the absolute value of the current I\ ) is 
maximized by working around the resonances e/ ~ fis = 



VII. DISCUSSION AND AHARONOV-BOHM 
OSCILLATIONS 

We have seen that there are two competing processes 
of currents, one where the two electrons proceed via dif- 
ferent dots into different leads, and one where the two 
electrons proceed via the same dot into the same lead. 
We will show now that these two processes also lead to 
different current oscillations in an Aharonov-Bohm loop 
which is threaded by an external magnetic flux (f>. For 
this let us consider now a setup where the two leads 1 
and 2 are connected such that they form an Aharonov- 
Bohm loop, where the electrons are injected from the 
left via the superconductor, traversing the upper (lead 
1) and lower (lead 2) arm of the loop before they rejoin 
to interfere and then exit into the same lead, where the 
current is then measured as a function of varying flux 
cf>. It is straightforward to analyze this setup with our re- 
sults obtained so far. In particular, each tunneling ampli- 
tude obtains a phase factor, Td 1 l 1 — * Tb 1 l 1 ^^ 2 ^ , and 
Td 2 l 2 -> T D2L2 e~^ /200 , where (j) = h/e is the single- 
electron flux quantum. For simplicity we also assume 
that the entire phase is acquired when the electron hops 
from the dot into the leads, so that the process dot-lead- 
dot gives basically the full Aharonov-Bohm phase factor 
e ±i<f>/<f> f the loop (and only a negligible amount of phase 
is picked up along the path from the superconductor to 
the dots). Now, we repeat the calculations of the transi- 
tion amplitude and find it to be of the following structure 
(/|T |») ~ T DlLl T D . 2L2 + T^ iLi e^o + T^e" 4 */*'. 
Here, the first term comes from the process via differ- 
ent leads (see (5.12)), where no Aharonov-Bohm phase is 
picked up. The Aharonov-Bohm phase appears in the 
remaining two terms, which come from processes via 
the same leads, either via lead 1 or lead 2 (see (6.5) 
and ( p-q )). The total current / is now obtained from 
|(/|T |z)| 2 , giving I = h + h+I AB , and the flux- 
dependent Aharonov-Bohm current Iab is given by 



I A B = V8hhF(e l ) cos 



F{ei 
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where, for simplicity, we have assumed that t\ = £2 = f-h 
and 71 = 72 = 7l- Here, the first term (different leads) is 
periodic in cj>g like for single-electron Aharonov-Bohm in- 
terference effects, while the second one (same leads) is pe- 
riodic in half the flux quantum </>o/2, describing thus the 
interference of two coherent electrons (similar single- and 
two-particle Aharonov-Bohm effects occur in the,Joseph- 
son current thr ough an Aharonov-Bohm loopi-l). It is 
clear from (7.1) that the h/e oscillation comes from the 



interference between a contribution where the two elec- 
trons travel through different arms with contributions 
where the two electrons travel through the same arm. 
Both Aharonov-Bohm oscillations with period h/e, and 
h/2e, vanish with decreasing I2, i.e. with increasing on- 
site repulsion U and/or gap A. However, their relative 
weight is given by yj I\/ 12, implying that the h/2e oscil- 
lations vanish faster than the h/e ones. This behavior 
is quite remarkable since it opens up the possibility to 
tune down the unwanted leakage process ~ I2 cos (2<f)/(f>o) 
where two electrons proceed via the same dot/lead by in- 
creasing U with a gate voltage applied to the dots. The 
dominant current contribution with period h/e comes 
then from the desired entangled electrons proceeding via 
different leads. On the other hand, if 1//1//2 < 1, which 
could become the case e.g. for kp5r > £/7> we are left 
with h/2e oscillations only. Note that dephasing pro- 
cesses which affect the orbital part suppress Iab- Still, 
the flux-independent current I\ + I2 can remain finite 
and contain electrons which are entangled in spin-space, 
provided that there is only negligible spin-orbit coupling 
so that the spin is still a good quantum number. 

We would like to mention another important feature of 
the Aharonov-Bohm effect under discussion, namely the 
relative phase shift between the amplitudes of tunneling 
to the same lead and to different leads, resulting in the 
addition al p rcfactor F(e{) in the first term of the rhs. of 
the Eq. (7.1). This phase shift is due to the fac t tha t there 



is a two-particle resonance in the amplitude (5-12) while 
t her e is onl y a single-particle resonance in the amplitudes 
( |6.5[ ) and ( |6.6| ) (we recall that the second resonance is 
suppressed by the Coulomb blockade effect). Thus, when 
the chemical potential us of the super condu ctor crosses 
the resonance, |e;| < 7l 



the amplitude ( p. 12 ) acquires an 



extra phase factor e 1 * 7 ', where <j) r — arg[l/(e; — 17^/2)]. 
Then the interference of the two amplitudes leads to the 
prefactor F(ei) = cos^v in the first term in the rhs of 
( |7.l|) . In particular, exactly at the middle of the reso- 
nance, ei = 0, the phase shift is <fr r = 7r/2, and thus the 
h/e oscillations vanish, since F(0) = cos(-7r/2) = 0. Note 
however, that although F = ±1 away from the resonance 
(|e; I ^> 7^) the h/e oscillations vanish again, now because 
the current I\ ~ ej s ^/L/ef vanishes. Thus the optimal 
regime for the observation of the Aharonov-Bohm effect 
is \ei\ ~7i. 

Finally, the preceding discussion shows that even if 
the spins of two electrons are entangled their associ- 
ated charge current does not reveal this spin-correlation 
in a simple Aharonov-Bohm interference experimcntE3. 
Only if we consider the current-current correlations 
(noise) in a beam splitter setup, can we detect also this 
spin-correlation in the transport current via its charge 
properticsE£. 



VIII. CONCLUSION 

We have proposed an entangler device that can cre- 
ate pairwise spin-entangled electrons and provide coher- 
ent injection by an Andreev process into different dots 
which are tunnel- coupled to leads. The unwanted pro- 
cess of both electrons tunneling into the same leads can 
be suppressed by increasing the Coulomb repulsion on 
the quantum dot. We have calculated the ratio of cur- 
rents of these two competing processes and shown that 
there exists a regime of experimental interest where the 
entangled current shows a resonance and assumes a finite 
value with both partners of the singlet being in different 
leads but having the same orbital energy. This entangler 
then satisfies the necessary requirements needed to detect 
the spin entanglement via transport and noise measure- 
ments. We also discussed the flux-dependent oscillations 
of the current in an Aharonov-Bohm loop. 

After finishing this prxyect we learned of a recent 
preprint by Lesovik et alE3 which also makes use of the 
Andreev process to generate spin-entangled electrons in 
the normal regime. There, the electrons are assumed to 
tunnel into fork-shaped normal leads with no Coulomb 
blockade behavior. The separation of the entangled elec- 
trons is achieved via energy filters so that the two elec- 
trons enter their corresponding lead with different orbital 
energies (while in the setup proposed here the energies 
can be made equal). 
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APPENDIX A: SUPPRESSION OF VIRTUAL 

STATES WITH BOTH ELECTRONS IN THE 

LEADS 

We have stated in the main text that the contributions 
of virtual states where two electrons are simultaneously in 
the leads are negligible. Here we estimate this contribu- 
tion and show that indeed it is suppressed by 7l/A^ <C 1 
(here the spin of the electrons is not important, and we 
set 71 = 72 = 7l for simplicity). First we consider the 
dominant transition from \DD) back to \DD) with the 
tunneling of only one electron to the lead, i.e. a sequence 
of the type \DD) -» \LD) -> \DD), and we find for the 



amplitude (cf. Eqs. (5.7,5.8)) 
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We compare this amplitude Adl with the amplitude 
All of the lowest-order process of tunneling of two 



electrons via the virtual state \LL), where both elec- 
trons are simultaneously in the leads, i.e. the sequence 
\DD) -v \LD) -> \LL) -> \DL) -v \DD). From now on 
we impose the resonance condition t\ = 0, and we find 
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where the first term in the bracket results from the se- 
quence of, say, electron f tunneling into lead f , then elec- 
tron 2 tunneling into lead 2, then electron 2 tunneling 
back into dot 2, and finally electron 1 tunneling back 
into dot f . While the second term in the bracket results 
from the sequence where the order of tunneling back to 
the dots is reversed, i.e. electron 1 tunnels back to its dot 
before electron 2 does. Note that due to this two terms 



in the bracket the two-particle pole in (A2) cancels 
Replacing £ k (. . .) with v L f^ A/x 



de(. 



we can write 
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Thus, comparing A^l with All, we see that indeed a 
virtual state involving two electrons simultaneously in 
the leads is suppressed by a factor of 7l/A/x compared 
to the one with only one electron in the leads. 



APPENDIX B: ELECTRON HOLE PAIR 
EXCITATION 

In this Appendix we consider a tunnel process where 
the two electrons starting from the superconductor tun- 
nel over different dots but during the process of repeated 
tunneling from the dots to the leads and back to the dots 
an electron from the Fermi sea hops on one of the dots 
(say dot 1) when this dot is empty. In principle, such con- 
tributions could destroy the desired entanglement since 
then a "wrong" spin can hop on the dot and the electron 
on the other dot (dot 2) would no longer be entangled 
with this electron (while the original partner electron dis- 
appears in the reservoir provided by the Fermi sea). We 
show now that in the regime A/i > 7; such electron-hole 
pair processes due to the Fermi sea are suppressed. We 
start with our consideration when the two electrons, after 
the Andreev process, are each on a different dot forming 
the I D-D) -state (we neglect spin in this consideration for 
simplicity). Instead of the amplitude (pq\T'\DD) calcu- 
lated in (|5.6|) we consider now the following process 
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The new sequence of interest in (Bl) is the amplitude 
containing the sum over n. For instance, let us consider 
the n = term, (DD\( i } H H DlLl ) 2 \DD), where we as- 
sume that the electron-hole excitation occurs in, say, lead 
1. From \DD), the tunnel Hamiltonian Hd 1 l 1 takes the 
electron from dot 1 to the state k in lead 1. Instead of 
tunneling back of this electron to dot 1 , an electron from 
the state k' with energy ek' < — A/i from the Fermi sea of 
lead 1 hops on dot 1. Now the dot-lead system is in the 



state \DD) = d^d^aik' a\^\i) . The sum over n resums the 
hoping back and forth of electron 2 from D2 to D2 , result- 
ing in the replacement of 77 in HutLi (Jti—HoI Hdi Li by 



7l/2. We perform the further resummation in (Bl) with 
this Fermi sea electron on dot 1 and the other electron 
still on dot 2, assuming that electron 1 in the state k in 
lead 1 is in its final state (and not a virtual state) . All the 
resummation processes in (Bl) are similar to the ones al- 



ready explained in the main text, except for having now 
an excitation with energy £k~ £k' > 0. The final state \pq) 
consists of two electrons in the lead states p and q (their 
multiple tunneling is resummed in T") and of the exci- 
tation with energy ek — £k', so \pq) = ai p a2q a ik' a ikN)- 
The normalized correction to the current, I e h/Ii, can be 
obtained by summing \A e h\ 2 /Ii over the final states \pq), 
and thus we arrive at the following integral for ei = 0, 
retaining only leading terms in jl/A/j,, and using energy 
conservation, ek — £k' + £p + £q =0, 
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We evaluate the integral in leading order and find 
- ^ 2^ { A-,) l0g {-) • 



(B2) 
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We see now that the current involving an electron-hole 
pair, I e h, is supp ressed compared to the main contribu- 
tion I\ (see (5.15)) by a factor of (-jl/A/i) 2 . 
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